What Is Claimed Is: 



1 LA method for bounding the solution set of a system of linear 

2 equations Ax = b, wherein A is an interval matrix and b is an interval vector, the 

3 method comprising: 

4 preconditioning the set of linear equations Ax = b by multiplying through 

5 by a matrix B to produce a preconditioned set of linear equations Mqx = r, 

6 wherein Mo = BA and r = Bb; 

H 7 widening the matrix Mo to produce a widened matrix M, wherein the 

p 8 midpoints of the elements of M form the identity matrix; and 

9 using M and r to compute the hull h of the system Mx = r, v^ich bounds 

* 1 0 the solution set of the system MqX = r. 

] 2, The method of claim 1 , wherein the method further comprises 

f^,^ 2 computing the matrix B by: 

H 3 computing an approximate center Ac of the matrix A; and 

^1 4 forming B by computing an approximate inverse of Ac, B = (Ac)" . 

1 3 . The method of claim 1 , wherein using M and r to compute the hull 

2 h involves: 

3 forming P as an inverse of the left endpoint of M; 

4 forming c^ = 1/(2P„ - 1) for z = 1, . . w; 

5 forming = (mfin)-\-sup(rj))Pu - e,^Vsup(rX wherein ei^ is a unit vector in 

6 which the z-th element is 1 and other elements are 0; 

7 setting inf{k) = c,z^ if > 0; 

8 setting inj\}i^ = if <. 0; and 

9 setting 5'w/?(h) = Ysupix). 
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1 4. The method of claim 1 , further comprising assuring that sup{r,) ^ 0 

2 by changing the sign of r, (and x,) if necessary. 

1 5 . The method of claim 1 , further comprising: 

2 determining if M is regular; and 

3 using the Gauss-Seidel process for computing the hull h if M is not 

4 regular. 

1 6. A computer-readable storage medium storing instructions that 

2 when executed by a computer cause the computer to perform a method for 

3 bounding the solution set of a system of linear equations Ax = b, wherein A is an 

4 interval matrix and b is an interval vector, the method comprising: 

5 preconditioning the set of linear equations Ax = b by multiplying through 

6 by a matrix B to produce a preconditioned set of linear equations Mqx = r, 

7 wherein Mo = BA and r = Bb; 

8 widening the matrix Mo to produce a widened matrix M, wherein the 

9 midpoints of the elements of M form the identity matrix; and 

1 0 using M and r to compute the hull h of the system Mx = r, which bounds 

1 1 the solution set of the system Mox = r. 

1 7. The computer-readable storage mediixm of claim 6, wherein the 

2 method further comprises computing the matrix B by: 

3 computing an approximate center Ac of the matrix A; and 

4 forming B by computing an approximate inverse of Ac, B = (Ac)'^ • 
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1 8. The computer-readable storage medium of claim 6, wherein using 

2 M and r to compute the hull h involves: 

3 forming P as an inverse of the left endpoint of M; 

4 forming = 1/(2?,, - 1) for ? = 1 , . . /a; 

5 forming zt = iinf{n)+supin))Pi^ - e^^?supir\ wherein ^/is a unit vector in 

6 which the z-th element is 1 and other elements are 0; 

7 setting infihj) = c^, if z^ > 0; 

8 setting inf{k) = if < 0; and 
H 9 setting sup{h) = Vsup{r), 

f] ] 9. The computer-readable storage medium of claim 6, wherein the 

=5 2 method fiirther comprises assuring that sup{n) > 0 by changing the sign of r, 

j; 3 (and Xj) if necessary. 



1 1 0. The computer-readable storage medium of claim 6, wherein the 

2 method further comprises: 

3 determining if M is regular; and 

4 using the Gauss-Seidel process for computing the hull h if M is not 

5 regular. 

1 1 1 . An apparatus that bounds the solution set of a system of linear 

2 equations Ax = b, wherein A is an interval matrix and b is an interval vector, 

3 comprising: 

4 a preconditioning mechanism that is configured to precondition the set of 

5 linear equations Ax = b by multiplying through by a matrix B to produce a 

6 preconditioned set of linear equations Mqx = r, wherein Mo == BA and r = Bb; 
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7 a widening mechanism that is configured to widen the matrix Mo to 

8 produce a widened matrix M , wherein the midpoints of the elements of M form 

9 the identity matrix; and 

10 a hull computing mechanism that is configured to use M and r to compute 

1 1 the hull h of the system Mx = r, which bounds the solution set of the system 

12 Mox = r. 

1 1 2. The apparatus of claim 1 1 , wherein the preconditioning mechanism 

2 is configured to: 

3 compute an approximate center Ac of the matrix A; and to 

4 form B by computing an approximate inverse of Ac, B = (Ac)'^ 

1 13. The apparatus of claim 1 1 , wherein the hull computing mechanism 

2 is configured to: 

3 form P as an inverse of the left endpoint of M; 

4 form - 1/(2P„ - 1) for z = 1, . . n; 

5 form z, = (mJ(r,)'\-sup(r>i)Pu - e^Ysup{r\ wherein e^i^ a unit vector in 

6 which the z-th element is 1 and other elements are 0; 

7 set infih) = c,z, if > 0; 

8 set infih^) = z^ifz^ < 0; and to 

9 set sup{h) = ¥sup(r), 

1 14. The apparatus of claim 1 1 , wherein the preconditioning mechanism 

2 is configured to assure that sup(r,) > 0 by changing the sign of (and xi) if 

3 necessary. 
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1 1 5 . The apparatus of claim 1 1 , wherein the preconditioning mechanism 

2 is configured to: 

3 determine if M is regular; and to 

4 terminate the process of computing the hull h if M is not regular. 

1 1 6. A method for bounding the solution set of a system of linear 

2 equations Ax = b by multiplying through by the matrix B to produce a 

3 preconditioned set of linear equations Mqx = r, wherein Mo = BA and r = Bb, the 
U 4 method comprising: 

5 assuring that supir;) > 0 by changing the sign of (and x,) if necessary; 

6 widening the matrix Mo to produce a widened matrix M, wherein the 
+ 7 midpoints of the elements of M form the identity matrix; and 

2 8 using M and r to compute the hull h of the system Mx = r, which bounds 

JU 9 the solution set of the system Mqx = r by, 

f* 10 forming P as an inverse of the left endpoint of M, 

P 11 forming c,= l/(2P„-l) for i=l, 

g 12 forming z, = iinf(r,)+sup(rJ)P„ - ejYsupir), wherein e.^is a 

13 unit vector in which the z-th element is 1 and other elements are 0, 

14 setting inj{h^) = c,z, if > 0, 

1 5 setting inf{h) ^ if z, < 0, and 

1 6 setting sup(yi) = Ysup{Y) . 

1 17. The method of claim 1 6, further comprising: 

2 determining if M is regular; and 

3 using the Gauss-Seidel process for computing the hull h if M is not 

4 regular. 
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1 18. The method of claim 1 6, wherein the method further comprises 

2 computing the matrix B by: 

3 computing an approximate center Ac of the matrix A; and 

4 forming B by computing an approximate inverse of Ac, B ^ (Ac)'^ 
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